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GUIDING  OF  ELECTRON  BEAMS  BT  HOLLOW  CONDUCTING  CHANNELS 


1.  Introduction 

An  electron  beam  propagating  inside  a  metallic  pipe  experiences 

magnetic  centering  forces  caused  by  return  currents  induced  in  the  pipe 

wall.  This  has  led  to  the  suggestion  that  an  annular  conducting  channel, 

created  by  a  laser  or  other  means,  could  guide  and  possibly  stabilize 

1—3 

relativistic  electron  beams  propagating  through  dense  gases.  A 
potential  application  is  in  inertial  confinement  fusion  reactors  where 
multiple  beams  from  isolated  diodes  propagate  through  several  meters  of  a 
dense  shielding  gas  toward  a  common  target.^  Advantages  of  hollow-channel 
guidance  over  other  proposed  schemes*-7  are  that  it  is  electrodeless  and 
requires  no  additional  power  supplies,  it  is  compatible  with  the  need  for  a 
dense  shielding  gas,  it  improves  beam  stability,  it  can  be  used 
repetitively,  and  the  laser  light  producing  the  channel  can  lie  outside  the 
target  and  need  not  interact  with  it. 

In  this  paper  we  provide  a  theoretical  basis  for  hollow-channel 
guidance  and  for  the  major  processes  affecting  it.  The  analysis  is 
restricted  to  gas  densities  sufficiently  high  that  a  scalar  conductivity 
adequately  describes  the  channel  electrical  properties.  To  understand  the 
interaction  between  a  beam  and  hollow  channel,  we  consider  a  rigid-rod  beam 
that  is  displaced  from  the  channel  axis  but  propagates  inside.  The 
displaced  beam  experiences  a  transverse  force  that  is  termed  "tracking"  if 
the  force  acts  to  center  the  beam  in  the  channel,  and  is  termed 
"detracking"  otherwise.  The  magnitude  and  duration  of  the  force  are 
computed  analytically  and  are  used  to  estimate  the  channel  conditions 
needed  for  successful  guidance.  Several  channel  configurations  (thin  and 
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thick  annuli,  rarefied  annuli,  and  a  magnetic  cage)  are  considered. 
Preliminary  computer  simulations  are  presented  to  support  the  analysis. 

All  of  these  tracking  schemes  are  passive  in  that  they  rely  on  channel 
currents  induced  by  the  beam.  Not  discussed  is  active  magnetic  guidance 
from  channel  currents  that  are  externally  driven.  The  dynamical  response 
of  the  beam  to  the  channel  forces  will  be  presented  in  a  later  paper. 

2.  Tracking  in  a  Metallic  Pipe 

To  illustrate  the  nature  and  magnitude  of  magnetic  tracking,  consider 

a  beam  propagating  inside  a  perfectly  conducting,  hollow  metallic  pipe. 

Injection  of  the  displaced  beam  into  the  pipe  induces  an  image  current  and 

charge  on  the  pipe.  The  pipe  return  current  I£  magnetically  repels  the 

beam  current  1^,  while  the  distributed  pipe  charge  Qc  electrostatically 

attracts  the  distributed  beam  charge  1^/00.  Here  0c  is'  the  axial  beam 

velocity  which  is  taken  to  be  everywhere  constant.  If  the  pipe  is 

evacuated  and  infinitely  long  with  no  conducting  end  plates,  Qc  ■*  lc/0c 

and  the  electrostatic  forces  overcome  the  magnetic  forces  to  order 
-2  2 

Y  =1-0  <<  1.  The  net  transverse  force  on  the  beam  is  then  weak  but 

detracking  and  propels  the  beam  into  the  pipe  wall.  Self-forces  similarly 
cause  the  beam  to  slowly  expand. 

If  the  pipe  is  filled  with  dense  gas,  the  beam  will  ionize  the  gas  and 
generate  conductivity  in  the  immediate  vicinity  of  the  beam.  This 
conductivity  shorts  out  the  electrostatic  fields  and  forces  on  the  beam, 
while  simultaneously  freezing  in  the  magnetic  fields  and  forces.  As  a 
result,  the  beam  begins  to  magnetically  pinch  and  track  the  pipe  center 
within  a  charge  neutralization  time  given  by 

T 

e 


l/4itab. 


(1) 


Beam  ionization  causes  to  rise  rapidly  so  that  the  charge-neutralization 
tine  for  an  intense  beam  in  a  dense  gas  is  typically  subnanosecond.  Ve 
shall  henceforth  generally  ignore  electrostatic  effects  unless  explicitly 
stated  otherwise. 

To  determine  the  magnitude  of  the  magnetic  tracking  force,  consider  a 
beam  of  fixed  radius  a^  that  lies  wholly  inside  a  pipe  of  inner  radius 
ac  >  a^.  The  beam  is  displaced  to  an  off-axis  location  y  <  a^a^.  In 
addition  to  the  pipe  current  Ic,  the  beam  induces  a  return  plasma  current 
Ip  in  the  ionized  gas.  If  the  beam  and  plasma  currents  are  assumed  to  be 
azimuthally  symmetric  about  the  displaced  axis  r^  -  y ,  the  axial  vector 
potential  from  these  currents  is  given  by 


a2i<El> 


In 


Er* 


(2) 


where  d  is  a  scale  length  determined  by  the  boundary  conditions  and 
I(  |r^-yj)is  the  net  current  flowing  within  a  radius  |r^-yj  about  the  y 
axis.  Observe  that  the  beam-centered  plasma  current  Ip  produces  no  net 
tracking  force  on  the  beam,  although  it  does  alter  the  magnetic  pinch  force 
on  the  beam.  Any  asymmetry  in  the  plasma  current  would,  of  course,  produce 
additional  tracking  or  detracking  forces.  We  show  in  Sec.  3  that  these 
forces  are  often  small  provided  the  plasma  conductivity  remains  centered 
about  the  beam.  Forces  from  the  plasma  currents  become  important,  however, 
once  the  beam  begins  to  track  (or  hose). 

Because  the  beam  and  plasma  currents  reside  off  axis,  the  induced  pipe 


current  I£  is  distributed  asymmetrically.  This  asymmetry  gives  rise  to  a 

magnetic  centering  force  on  the  beam  that  can  be  computed  by  replacing  the 

pipe  by  an  image  line-current  Ic  located  at  r^  =  y  .  The  image  line- 
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current  generates  an  axial  vector  potential  given  by 


(3) 


A*2(£l> 


where  d'  is  another  normalization  length.  The  pipe  boundary  condition  for 
the  total  vector  potential  is 


zl 


+  Az2  "  0 


(A) 


at  r  =  ac«  Combining  Eqs.  (2)-(4)  produces,  after  some  algebraic 
manipulation,  the  following  results: 


I  =  - 
c 


(5) 


X'  =  <ac/y)2y  , 


(6) 


and 


d'/d  =  ac/y 


(7) 


where  I  =  I,  +  I  is  the  net  current.  The  total  current,  I.  +  I  +  I  , 
n  b  p  -  b  p  c 

equals  zero. 

The  image  line-current  Ic  located  at  y'  repels  the  beam  current  1^ 
located  at  y  with  a  distributed  force  given  by 


Nh  Fm  .  K  -4-S. 

b  ~m  2  i  . 

c  k -i 


,  ,  .  2  Ib1n 

2  X)  =  "  2  2  2  X* 

I  c  ac-y 


(8) 


The  average  magnetic  tracking  force  per  beam  electron  is  thus 


2ef5  n 

=  c  ~2 — 2  X 
a  -y 
c  J 


F 


(9) 


where  *  -  I^/ePc  is  the  beam  line  density.  Result  (9)  characterizes  the 
maximum  magnetic  tracking  force  that  a  hollow  conducting  channel  can  exert 
on  a  beam.  Note  that  Fm  attains  a  magnitude  comparable  to  the  average 
monopole  pinch  force,  Fq  =  efH^/a^c,  as  the  displacement  y  ■*  ac-a^. 

Although  the  transverse  beam  force  is  tracking,  the  transverse  force 
on  the  counter-flowing  plasma  current  1^  is  detracking.  In  a  dense  gas  the 
plasma  channel  does  not  physically  displace,  but  the  plasma  electrons 
instead  establish  a  Hall  electrostatic  field  to  offset  the  magnetic 
detracking  force: 


2v,  1 

p _ 1  _ n_ 

~H  2  2  2  ^ 

c  a„-y 


(10) 


where  v^  is  the  axial  plasma-electron  drift  velocity.  This  Hall  field 
produces  a  weak  detracking  force  on  the  beam  electrons  which  can  be  ignored 
to  order 


eEH/Fm  =  vd/0c  «  1.  (11) 

The  force  would  be  tracking  if  the  beam  and  plasma  current-carriers  were 
oppositely  charged.  The  Hall  effect  also  produces  a  slight  lateral  shift 
in  the  plasma  current  that  has  a  yet  weaker  influence  on  tracking. 

3 .  Tracking  in  a  Conducting  Gaseous  Annulus 

Replacing  the  metallic  pipe  with  a  preionized  and  conducting  gaseous 
annulus  introduces  three  new  features.  Two  arise  from  the  finite 


conductance  of  the  annulus.  The  third  arises  from  the  tenuousness  of  the 
annulus  which  permits  beam  electrons  to  easily  penetrate  it. 


The  first  new  feature  is  that  electrostatic  fields  between  the  inner 


and  ouer  radii  of  the  annulus  are  no  longer  fully  suppressed.  These 

fields  reduce  the  return  current  carried  by  the  annulus  and  thereby  reduce 

the  magnetic  tracking  force  on  the  beam.  Although  the  channel  charge  and 

corresponding  electrostatic  detracking  force  are  similarly  reduced,  an 

additional  detracking  force  arises  from  bipolar  charging  of  the  channel. 

Positive  and  negative  charge  sheaths  form  on  the  inner  and  outer  channel 

edges,  respectively,  to  produce  an  electrostatic  detracking  force  that  is 

o 

the  analogue  of  the  tracking  force  produced  by  solid  channels.  A  very  low 
conductivity  in  a  hollow  channel  thus  produces  a  net  beam  force  that  is 
initially  detracking  although  weak. 

The  second  new  feature  is  that  magnetic  flux  diffuses  through  the 

resistive  annulus.  Any  azimuthal  asymmetry  in  the  channel  current  thus 

relaxes  on  a  finite  channel  dipole  decay  time  Tc^.  Because  a  symmetric 

channel  current  produces  no  fields  or  forces  on  beam  electrons  inside  the 

channel,  the  channel  tracking  force  relaxes  on  Tcl>  The  integrated  channel 

current  Ic  relaxes  on  a  channel  monopole  time  x^q  which  is  usually  longer, 

x  >  x  . . 
co  cl 

Spreading  of  the  beam  through  the  gaseous  channel  is  the  third  new 
feature.  A  circular  ring  of  beam  current  which  lies  entirely  outside  the 
channel  produces  zero  electromagnetic  field  inside  the  ring,  due  to  the 
assumed  azimuthal  symmetry  of  the  beam.  Such  a  ring  thus  experiences  zero 
net  force  from  the  channel.  A  ring  which  passes  through  the  channel  always 
experiences,  however,  a  detracking  force  caused  by  channel  return  current 
flowing  outside  the  ring;  see  Fig.  1.  The  net  magnetic  force  on  a  beam 
with  an  extended  radial  profile  is  thus  comprised  of  a  tracking  component 
from  beam  rings  inside  the  channel,  and  a  detracking  component  from  beam 
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rings  intersecting  the  channel.  At  late  times  only  the  detracking 
component  survives.  The  detracking  effect  from  the  beam  wings  was  ignored 
in  deriving  Eq.  (9). 

These  new  features  suggest  that  a  conducting  annular  channel  can 
strongly  guide  an  electron  beam  only  if  three  conditions  are  met.  First, 
the  channel  should  have  sufficient  conductivity  a to  suppress  internal 
electrostatic  fields:  typically, 

4n<x  a  /c  >  1  (12) 

c  c  v  ' 

where  ac  is  the  inner  channel  radius.  Otherwise,  electrostatic  detracking 
of  the  beam  head  may  take  place. 

Second,  the  magnetic  tracking  force  from  the  channel  current  should 
persist  longer  than  the  electrostatic  detracking  force  from  the  channel 
charge.  Strong  magnetic  tracking  can  thus  occur  only  if  the  channel 
magnetic  dipole  decay  time  exceeds  the  beam  charge-neutralization  time, 


cl 


t  . 
e 


At 


the  beam  pinch  point  where 
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(13) 


4nababp/c  = 


(14) 


condition  (13)  reduces  with  the  help  of  definition  (1)  to 


Tcl  >  abp/c- 


(15) 


Here  a^p  is  the  beam  radius  at  the  pinch  point. 

The  third  condition  is  that  the  annulus  should  lie  outside  the  beam  to 


minimize  magnetic  detracking  by  the  beam  wings: 


If  the  latter  condition  is  not  met,  the  hollow  channel  can  act  as  a  solid 
channel  and  eject  the  beam  via  magnetic  repulsion  from  the  opposing  channel 
current.  Increasing  the  inner  channel  radius  to  much  beyond  requirement 
(16)  decreases,  however,  the  tracking  force  as  Eq.  (9)  demonstrates. 

Optimal  beam  guidance  is  therefore  expected  to  occur  when 

1  <  a  /a.  <2.  (17) 

c  bp 

When  conditions  (12),  (15)  and  (16)  are  satisfied,  the  net  force  on 

the  beam  head  is  tracking  and  remains  tracking  (in  dense  gases)  to  beyond 

the  pinch  point.  After  the  pinch  point,  magnetic  coupling  forces  become 

strong  and  can  overcome  any  residual  detracking  force  arising  from  beam- 

channel  overlap.  The  coupling  forces  thus  cause  the  beam  body  to  follow 

the  head.  These  same  coupling  forces  also  give  rise  to  the  resistive  hose 
11  12 

instability  ’  briefly  discussed  in  Sec.  3.4. 


3.1  Linearized  Analysis 

Investigating  the  consequences  of  condition  (15)  requires  computation 
of  the  channel  dipole  decay  time  To  compute  this  time  and  to  analyze 
the  problem  in  further  detail,  we  use  Ampere's  law  for  the  axial  vector 


potential  A : 


V^A  = 
1  z 


,  ,  9A  . 

c“  (Jb  "  IT) 


(18) 


where  is  the  beam  current  density,  a  is  the  plasma  conductivity,  and 
t  a  ct-z  is  the  distance  behind  the  beam  head.  Here  we  have  employed  both 
the  ul t rarelat i vis t ic  frozen  approximation, 


(19) 


and  the  magnetostatic  approximation, 

*  3A, 

Ez  -  -  h  (Az " *>  “ "  w  (20) 

where  E  is  the  axial  electric  field  and  ♦  is  the  electrostatic  potential, 
z 

The  neglect  of  ♦  is  strictly  justified  only  within  the  annulus  itself, 
where  4ne  a  /c  >  1,  and  behind  the  pinch  point,  where  4neha./c  >  1.  Ahead 
of  the  pinch  point,  our  results  for  the  magnitude  and  duration  of  the 
magnetic  channel  tracking  force  still  apply,  at  least  approximately,  when 
expressed  in  terms  of  the  channel  current  Ic;  in  this  region,  however,  Ic 
depends  strongly  on  Not  computed  in  this  analysis  are  electrostatic 
effects  (such  as  detracking  by  channel  charge)  which  are  complicated  by  the 
dynamics  of  space-charge  neutralization. 

To  simplify  the  problem,  we  consider  small  displacements  y  «  a^  <  a£ 
and  divide  the  plasma  conductivity  into  two  components.  The  component  ob 
is  azimuthally  symmetric  about  the  beam  axis  r^  =  whH®  the  component  oc 
is  azimuthally  symmetric  about  the  channel  axis  r^  -  0.  As  discussed  in 
Secs.  4  and  5,  processes  such  as  avalanching  and  high-order  chemistry 
effects  can  complicate  and  modify  this  separation  of  a.  Here  we  ignore 
these  complications  and  perform  a  straightforward  linearized  azimuthal 
expansion  about  the  channel  axis: 

Az(r,0,O  =  A0(rt£)  +  Ax(r*0  cos9  +  ...  (21a) 

3J. 

Jb(r,e,0  =  Jb(r',C)  -  y  3^7  cos©  +  ...  (21b) 
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<jb(r,e,t)  =  <Tb(r',t)  -  y  jp-  cose  +  ... 


(21c) 


where  r'  =  |r^  -  jr|.  We  treat  the  dipole  terms  as  small  perturbations  and 
set  r'  -  r.  Inserting  expansion  (21)  into  Eq.  (18)  then  yields 


-  —  r  —  -  -  —  lJw  -  +  °c)  W 


r  3r  3r 


for  the  monopole  potential  Aq,  and 


„  .  3J.  3  a,  3A  3A. 

If  F  h  '*1  *  r  »  ST  -  *  5T  5T  *  <"b  *  «e>  ST 


for  the  dipole  potential  A^. 

Each  beam  electron  experiences  a  transverse  magnetic  force  given  by 


ee  ?iV 


The  average  tracking  force  for  the  entire  beam  is  thus 


l>:5 
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Fm  -  -  4J‘\  r; 


J.  3A  .  0  9 A 

-  .  b  ~  z  sine  z 

-  ep  dr  r  d9  cos9  jp  -  -p~  55- 
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where  y  is  the  unit  displacement  vector.  Substituting  expansion  (21)  into 
Eq.  (25)  produces 


dr  r  Jb(r  If  rAl) 


3J,  3A„ 
_ b  _ o 

3r  3r 


Integrating  both  terms  by  parts  and  using  the  dipole  boundary  condition, 


Li“  (J  9r  rAl)  "  °» 


r-x» 


(27) 


we  obtain 


!*[4b(lr£ 


(28) 


Here  the  beam  radial  profile  is  defined  by 


Ib<r»« 


dr  2 nr  Jb(r,C)/Ib(C). 


(29) 


Substituting  expressions  (22)  and  (23)  into  Eq.  (28)  produces 


2neg 


drjl 


3Jb  3%  3Ao  (  \  3A1 

y  dT  “  y  IT  dt~  +l  %  +  acJ  IT 


♦  y 


fit 

9r 


dA 

Jb  "  (ab  +  ®c)  W 


Integrating  the  first  two  terms  by  parts  reduces  this  to 


2ne6 
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dr 


dA,  dIL  dA. 


dA 


Vc  dT 
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y  F  ffc  aT  +  Vb  IF 


(30) 


(31) 


where 


dA 

Aib  ■  Ai  *  y  5r 


(32) 


Result  (31)  for  the  average  magnetic  deflection  force  Fm  has  the 

following  physical  basis.  The  first  term  in  the  integrand,  1^  ac  9Aj/9£, 

represents  the  contribution  from  channel  dipole  currents  that  produce  the 

channel  tracking  force  previously  described.  The  seco'nd  term, 

-y  91^/ 3r  ac  9Aq/9C,  represents  the  channel  detracking  force  caused  by 

overlap  of  the  beam  current  (J^  *  91^/ 3r)  and  the  channel  current 

(J  =  -  a  9A  /3C).  The  last  term,  I.  a.  3A,./9C,  arises  from  plasma 
cco  b  b  lb 

current  in  the  beam-centered  conductivity  o^.  This  plasma  current  produces 
a  net  transverse  force  on  the  beam  only  if  it  is  asymmetrically  distributed 
with  respect  to  the  beam.  Here  Alb  is  the  dipole  potential  for  an 
expansion  about  the  beam  axis  r^  =  jr.  This  potential  is  defined  by 
Eq.  (32)  and  satisfies 


9_  1  9_ 

9r  r  9r  rAlb 


Art 

c 


3ac  3Ao  ,  .  3Albl 

y  W  JC  +  (ffb  +  ffc>  9C~ 


(33) 


3.2  Thin  Annulus  ("Resistive  Pipe") 

The  magnitude  and  duration  of  the  tracking  and  detracking  components 

of  Fm  can  be  computed  using  Eqs.  (22),  (23),  and  (33).  To  further  simplify 

the  problem,  consider  first  a  narrow  channel  with  a  specified  distributed 

resistance  R  (O: 

c 

ac ( 1 • O  =  &(r  -  ac)/2nacRc(C)  (34) 

where  6  is  the  Dirac  delta  function.  Inserting  this  form  into  Eq.  (23), 
applying  boundary  condition  (27)  and  definition  (32),  and  performing  two 
radial  integrations  yields  for  the  channel  dipole  potential 
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w 


a  ® 

*  c  p 

m  — 

dr  r  dr 

« 

J  acc 

c 

r 

+  % 

3  (.  3Ao) 

3t  iAl  +  y  3r  J 

3A-  a  /  3A  v 

ffc  w +  y  If  K  -  %  ar) 


1  3  y  21<v<> 

-  rr  r  Ai(v;> + 1 - ^ — 


a 

r  C 


4lt 

acc 


dr  r 


dr  ffb  3C  Alb‘ 


(35) 


Here  the  net  current  flowing  just  inside  the  channel  is  given  by 


i<vo 


dr  2nr  I  J,  - 


(36) 


Inspection  of  the  term  on  the  left  and  the  first  term  on  the  right  of 
Eq.  (35)  indicates  that  Aj(ac>0  relaxes  on  a  characteristic  dipole  decay 
length  given  by 


CT  .  =  1/R  c. 
cl  c 


(37) 


The  channel  dipole  current  and  tracking  force  relax  similarly. 

Applying  the  same  technique  to  Eq.  (22)  yields  for  the  channel 
monopole  potential 


Vac’C) 


2  ln(b/ac) 

iTc 

c 


tc  V*c'°  *Hr I<r-°  <38> 
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where  b  >  a£  is  the  monopole  boundary  radius  and  I(r,C)  is  the  net  current, 
exclusive  of  the  channel  current,  flowing  within  radius  r: 

r 

'  3A 

I(r,C)  ■  dr  2 nr  -  crb  (39) 

o 

Inspection  of  Eq.  (38)  indicates  that  the  channel  monopole  potential, 
field,  and  current  relax  on  a  characteristic  monopole  decay  length  given  by 

ctco  =  (2/Rcc)ln(b/ac).  (40) 

The  channel  monopole  and  dipole  decay  times  are  thus  related  by 

t  / t  .  =2  ln(b/a  ).  (41) 

CO  cl  c 

The  boundary  radius  b  in  these  equations  is  where  the  monopole 
electric  field  falls  to  zero: 

Ez0(b,O  =  0.  (42a) 

In  the  magnetostatic  limit  considered  here,  we  impose  the  equivalent 
condition  that 

Aq  =  0  (42b) 

for  all  r  >  b.  In  the  absence  of  external  return-current  structures,  as 
was  assumed  for  the  dipole  boundary  condition  (28),  the  monopole  boundary 
radius  b  equals  the  vacuum  radius  b  where  the  conductivity  becomes  too 
low  to  neutralize  the  space-charge  fields.  Typically, 

4rtar/c  <0.1  (43) 

for  all  r  >  b  >  a  .  If  a  metallic  pipe  or  similar  structure  is  present 
-  vac  c 

at  a  radius  b  .  ,  the  boundary  b  equals  the  lesser  of  b  . _  and  b  _  . 

pipe  pip®  vac 

More  importantly,  the  dipole  boundary  condition  (28)  must  then  be  replaced 


with  the  requirement  that  Ai^pipe*^)  *=  0.  This  requirement  gives  rise  to 
an  additional  magnetic  tracking  force  caused  by  return  currents  induced  in 
the  pipe. 

Ve  have  thus  far  shown  that  t^e  channel  tracking  force,  as  represented 

by  the  first  term  in  the  integrand  of  Eq.  (31),  relaxes  on  a  channel  dipole 
2 

time  tcj  =  1/Rcc  .  The  second  term,  which  represents  a  detracking  force 
from  beam-channel  overlap,  relaxes  more  slowly  on  a  channel  monopole  time 
Tco  *  Tcl’  The  detracking  force  is  small,  however,  if  the  overlap  is 
small.  From  Eqs.  (31),  (35)  and  (38)  one  can  show  that  the  ratio  of  the 
overlap  detracking  force  to  the  tracking  force  is  initially  given  by 

lFco/Ftl  =  "VVV^VV^  1  1  (44) 

where  I^a^t)  is  the  beam  current  within  r  *  ac-  This  result  can  also  be 
inferred  by  using  Fig.  1  to  estimate  the  detracking  force  from  beam-channel 
overlap  and  comparing  it  to  the  tracking  force  given  by  Eq.  (9).  For  the 
infinitesimally  thin  channels  considered  here,  the  net  channel  tracking 
force  is  maximized  for  Bennett  or  flat-topped  beams  when  a£  -  a^.  For 
channels  of  finite  breadth,  maximum  tracking  is  generally  realized  when 
condition  (17)  is  satisfied.  If  ac  >  2a^p  +  y,  the  detracking  force  can 
usually  be  neglected  altogether. 

We  have  yet  to  assess  rV'  last  term  in  Eq.  (31)  which  represents  the 
dipole  force  arising  from  plasma  current  induced  in  the  beam-generated 
conductivity  o^.  This  component  of  the  tracking  force  was  neglected  in  the 
nonlinear  derivation  given  in  Sec.  2.  We  show  here  that  such  neglect  is 
often  well  justified,  provided  the  beam  is  rigidly  displaced.  This 
component  becomes  important,  however,  once  the  beam  begins  to  track  or 
hose. 
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■  ■  t 


is  the  beam-centered  resistance  per  unit  length.  Equations  (46)  and  (47b) 


indicate  that  relaxes  on  a  characteristic  dipole  decay  length  given  by 

2 

ct^j  =  1/R^c  =  ^^b^2*  (49) 

The  detracking  force  from  dipole  currents  in  the  beam-centered  conductivity 
)i.e.,  the  last  term  in  Eq.  (32))  thus  relaxes  on  a  time  “ 
(Rc/R^)tc^.  This  detracking  is  short-lived  and  can  be  neglected  if 
R^  »  Rc>  as  is  usually  true  in  the  beam  head. 

In  evaluating  the  characteristic  relaxation  times  for  and 
Aj(ac,t)>  ve  neglected  the  coupling  betveen  them.  Such  coupling  alters  the 
actual  relaxation  times  and  thereby  alters  the  tracking  force.  The  maximum 
effect  that  A^  has  on  A^(ac,C)  can  be  estimated  by  using  Eqs.  (46)  and 
(47b)  to  compute  an  approximate  upper  bound  for  3/3C  A^: 


f3Alb] 

(r/ac) 

N 

3A1(ac,C) 

ac 

/ 

~  "  2  2 
max  1-r  /2a. 

iRcj 

ac 

(50) 


Substituting  this  upper  bound  into  Eq.  (35)  produces 


1  3  y  2I<VC> 

ll(acC)  =  R  c  3C  Al(ac,C)  *  a  c 


4n 

+  2 
a  c 
c 


'V 

R 

c 


3C  Al(ac’C) 


dr  r 


dr  a 


b  ,  2/0  2 

!-r  /2ab 


(51) 


where  the  limits  in  the  integrands  have  been  changed  to  reflect  the 
assumption  that  -»  0  for  r  >  a^.  For  flat  we  can  readily  evaluate 
the  integrals  and  rewrite  Eq.  (51)  as 


2, 


1  -  Ot  —jr 


j  3A 1 ( ac  •  C) 


R  c 

c 


ac 


2I(ac,Q 

a  c 

c 


A1 (ac  * 


(52) 
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V.  A. 


where  a  *  ln(16)  -  2  =  0.77.  In  this  example,  A^(a  ,£)  thus  relaxes  on  a 

'  2 
time  t  ,  that  is  somewhat  shorter  than  T  .  *  1/R  c  : 
cl  cl  c 

Tcl/Tcl  -  1  -  a  (ab/ac)2>  °-23*  (53) 

9 

For  a  >  1.5  a,  ,  t  ,  >  T  ,/2.  Moreover,  if  R.  >>  R  as  is  usually  true  in 
c  —  b  cl  cl  be 

9 

the  beam  head,  x  ^  -»  t  ^  regardless  of  the  ratio  a^/a^- 

The  accelerated  relaxation  of  A^(ac,C)  increases  the  maximum  channel 
tracking  force  while  decreasing  its  duration.  The  increase  in  magnitude 
represents  the  reactive  response  of  the  conducting  channel  to  the 
detracking  current  induced  in  o^.  The  net  tracking  force  is  nonetheless 
weakened.  In  the  absence  of  beam-channel  overlap,  the  maximum  net  tracking 
force  is  given  from  Eqs.  (31),  (50),  and  (52)  by 


.(max) 


0.23 


1-0. 77a?/a2 

D  C 


2e|3yI(ac,C) 

2 

a  c 
c 


2,2, 


(54) 


Here  we  have  assumed  flat  (i.e.,  »  r  /a^)  for  r  <  a^.  Except  for  the 

term  in  square  brackets,  Eq.  (54)  represents  the  linearized  equivalent  of 
Eq.  (9)  in  the  limit  y  <<  ac>  The  term  in  brackets  represents  the  maximum 
reduction  in  the  tracking  force  from  dipole  currents  in  a^.  If  R^  >>  Rc, 
no  reduction  occurs,  and  the  term  in  brackets  should  be  set  to  unity. 


3.3  Thick  Annulus 

Consider  now  the  opposite  limit  of  a  broad,  hollow  conducting  annulus 

for  which  a  =  0  for  r  <  a  .  If  we  neglect  both  beam-channel  overlap 
c  c 

(a.  <  a  )  and  all  effects  from  dipole  currents  in  a,  ,  Eqs.  (23)  and  (31) 

DC  D 

reduce  respectively  to 


1, 


3  13  . 

dr  r  3r  rAl‘ 


3A 


b  3C 


(55) 


and 


F 

m 


2neg 

c 


I 


dr 


a 


(56) 


Here  ve  have  used  definition  (32)  to  rewrite  and  eliminate  all  terms 
containing  3Aj^/3C  from  Eq.  (23),  and  have  set  1^  -  1  for  r  >  ac  in  Eq. 
(31). 

By  integrating  Eq.  (55)  twice,  we  obtain  with  the  help  of  Eqs.  (27) 
and  (39)  for  all  r  <  ac, 


rA 


1 


3A  . 

5T) 


+ 


a 


3^ 

c  3C  J 


I(r,C) 


2  nr 


3A. 


J  "c  5T 


(57) 


where  we  have  used  the  fact  that  =  0  for  r  <  ac.  Combining  this  with 
Eq.  (56)  yields 


F  (C)  _  2epyl(rt(;) 
2 


e@A1(r,C) 

r 


(58) 


for  all  r  <  a^.  This  equation  demonstrates  that  diffusion  of  the  internal 

dipole  potential  A^  controls  the  evolution  of  the  channel  tracking  force. 

At  late  times,  A,  -»  2yI(r,C)/rc  while  F  -»  0. 

1  m 

To  estimate  the  characteristic  relaxation  rate  for  F  ,  we  take  a  to  be 

m  c 


p 

I# 


0  (r)  =  0  H(r-a  ) 
c'  '  co  v  c' 


(59) 


where  H  is  the  Heaviside  step  function.  The  solution  to  Eq.  (55)  for 
r  >  ac  and  for  constant  or  slowly  varying  I(ac,C)  is  then  given  by 


Ax(r,0 


2yI(ac,C)  -fl0cor2/cC 
rc  6 


(60) 


Because  must  be  everywhere  continuous,  Eqs.  (58)  and  (60)  may  be 
combined  at  r  =  a£  to  produce 


F_(  C) 


2e0yI(ac,C) 


2 

a  c 
c 


/  2  \ 
-no  a  /cC 


1-e 


co  c 


(61) 


The  characteristic  relaxation  length  for  the  tracking  force  produced  by  a 
broad  hollow  channel  of  inner  radius  a£  and  conductivity  ©co  is  thus 


CT 


cl 


0  /c. 

co 


(62) 


Note  that  the  tracking  force  in  broad  channels  does  not  follow  the  simple 
exponential  decay  [exp(-C/cTc^) ]  predicted  for  narrow  channels.  The  slow 
decay  in  broad  channels  occurs  because  the  dipole  channel  tracking  current 
diffuses  radially  outward  before  it  decays.  The  diffusion  is  reflected  in 
Eqs.  (60)  and  (61). 
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3.4  Optimal  Channel  Parameters 


3 


$ 


I 


a 


We  now  restate  the  channel  parameters  needed  to  effect  strong  magnetic 


tracking  at  the  beam  pinch  point.  We  consider  both  broad  and  narrow 


conducting  annuli  and  arbitrarily  assume  that  the  beam  radius  at  the  pinch 


point,  as  defined  by  condition  (14),  is  =  3  cm.  This  is  the  main  beam 


parameter  that  influences  the  choice  of  optimal  channel  parameters. 


The  optimal  inner  channel  radius  is  given  by  condition  (17),  which  is  a 


compromise  between  channel  tracking  and  the  detracking  caused  by  beam- 


channel  overlap.  A  reasonable  choice  is  a  =  5  cm  for  a.  =3  cm. 

c  bp 


The  minimum  channel  conductivity  needed  to  subdue  internal  channel 


8  -1 

electrostatics  is  given  by  condition  (12):  acQ  >  c/4rtac  =  5  x  10  s  for 


ac  =  5  cm.  This  typically  corresponds  to  a  channel  ionization  fraction  of 


n  /n  >  2  x  10  where  n  is  the  channel  electron  density  and  n  is  the 
ec  g  ~  ec  g 


gas  density.  The  gas  is  thus  weakly  ionized.  The  low  degree  of  ionization 


suggests  that  the  channel  could  be  created  using  a  laser  to  ionize  a  seed 


gas  added  to  an  inert  base  gas.' 


For  thin  annuli  of  width  &a c  <<  ac»  the  magnetic  conditions  (15)  and 


(37)  impose  an  upper  bound  on  the  channel  resistance: 


R  <  1/a.  c. 
c  bp 


This  condition  insures  that  magnetic  tracking  outlasts  and  dominates 


electrostatic  detracking.  For  a^  =  3  cm,  a  channel  resistance  of 


-11  2 

R  <10  s/cm  =  10  S/cm  is  sufficient.  To  obtain  low  R  for  &a  <<  a 
c  c  c  c 


requires,  however,  a  channel  conductivity  much  higher  than  stipulated  by 


condition  (12): 


(64) 


4ji a  a  /c  =  2/R  Sa  c  >  2a,  /Sa  >>  1. 
co  c  c  c  bp  c 


Such  high  conductivity  is  difficult  to  establish  and  maintain  in  dense 


gases.  Narrow  conducting  hollow  channels  are  thus  often  impractical  for 


tracking  (particularly  if  formed  by  lasers  or  similar  gas-ionization 


methods) . 


For  broad  channels  of  width  Sac  >  ac>  condition  (64)  is  replaced  by 


4 no  a  /c  >  4  a.  /a 

co  c  bp  c 


which  is  derived  using  Eqs.  (15)  and  (62).  In  the  present  example, 


condition  (65)  is  only  slightly  more  stringent  than  condition  (12)  and  is 
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satisfied  for  o  >  10  s  The  total  channel  resistance  still  satisfies 
co 


condition  (63).  Note  that  channel  conductivity  generated  outside  r  -  2a 


is  of  little  benefit  to  tracking  (although  it  can  aid  beam  stability). 


Satisfying  the  magnetic  condition  (15)  provides  two  additional 


benefits.  One  is  that  the  channel  dipole  decay  time  is  larger  than 


four  times  the  beam  dipole  decay  time  at  the  pinch  point 


t  ,  >  a./c  =  4t,  .  ( C  ) . 
cl  bp  blv  p' 


Magnetic  detracking  from  dipole  currents  in  can  therefore  be  safely 


neglected  up  to  and  somewhat  beyond  the  pinch  point. 


A  second  and  potentially  more  important  benefit  is  the  impact  that 


condition  (66)  has  on  beam  stability.  The  most  virulent  instability  is 


beam  hose  which  develops  when  a  nonrigid  beam  displacement  gets  out  of 


phase  with  the  induced  restoring  forces.  The  phase  delay  and  growth  rate 


for  this  instability  are  governed  by  the  dipole  decay  time.  A  hollow 


conducting  channel  increases  the  dipole  decay  time  and  thus  imposes  an 


upper  bound  on  the  hose  growth  rate.  In  the  beam  head  where  x  ^  >> 


hose  develops  as  a  resistive  vail  instability  with  a  growth  rate 

proportional  to  .  In  the  beam  body  where  becomes  large  and 

11  12 

>>  tcj,  the  usual  resistive  hose  instability  ’  develops  with  a 
growth  rate  proportional  to  .  A  hollow  channel  can  thus  reduce  hose 

growth  in  the  beam  head,  and  thereby  indirectly  reduce  hose  growth  (but  not 


the  growth  rate)  in  the  beam  body.  A  hollow  channel  additionally  lessens 
hose  instability  (and  axisymmetric  instabilities)  by  reducing  the 
destabilizing  plasma  return  current  in  the  immediate  vicinity  of  the  beam. 
Although  the  present  analysis  is  restricted  to  rigid  beams,  we  note  that 


the  stabilizing  influence  of  hollow  channels  has  been  observed  both 


4  n  1.13  .  .  ,  n  2,14 

numerically  and  experimentally. 


4 .  Alternative  Schemes 

14 

A  variation  on  the  resistive  pipe  is  a  magnetic  cage  consisting  of  a 
ring  of  three  or  more  solid  conducting  channels  encircling  the  beam.  A 
four-conductor  cage  is  depicted  in  Fig.  2.  The  azimuthal  asymmetry  of  the 
cage  strengthens  the  tracking  force  along  some  directions  but  weakens  it 
along  others.  At  late  times  when  the  channel  current  has  become  equally 
distributed,  the  force  becomes  detracking  for  beam  deflections  between 
adjacent  channels,  even  in  the  absence  of  beam-channel  overlap.  The  beam 
head  nonetheless  tracks  and  guides  the  beam  body.  The  cage  is  thus  an 
acceptable  alternative  to  smooth  annular  channels,  although  it  does 
slightly  degrade  tracking  and  requires  a  higher  degree  of  local  gas 
ionization  to  obtain  a  given  net  resistance  R^.  The  ability  of  a  cage  to 
guide  and  stabilize  a  beam  was  demonstrated  experimentally  by  Raleigh  and 


Fernsler  using  four  resistive  rods. 


•v\.  v ...  ... 


Creating  high  channel  ionization  for  either  the  pipe  or  cage  mode  can 


be  a  formidable  problem.  One  means  of  circumventing  this  problem  is  to 

2 

utilize  avalanching  in  the  space-charge  fields  of  the  beam.  The  space- 
charge  fields  peak  off  axis  and  are  capable  of  creating  high  conductivity 
outside  an  intense  beam.  If  a  rarefied  or  weakly  preionized  hollow  channel 
surrounds  the  beam,  preferential  avalanching  can  transform  the  channel  into 
a  good  conductor. 

The  degree  of  rarefaction  or  preionization  required  to  create  a 
conducting  channel  suitable  for  tracking  can  be  estimated  as  follows. 
Avalanching  causes  the  conductivity  outside  the  beam  to  grow  according  to 


c 


p  v  a 


(67) 


where  p  is  the  local  gas  density  and  v(E/p)  is  the  reduced  avalanche  rate. 
Here  we  have  dropped  the  convective  terms  (as  is  often  done  in  simulations) 
because  they  little  affect  the  final  conductivity.  Also  dropped  is  beam 
direct-production  on  the  assumption  that  •*  0  for  r  >  ac  >  a^. 

In  the  low-conductivity  regions  (4fiur/c  <<  1)  outside  the  beam,  the 
dominant  electric  field  is  the  radial  field  which  in  axisymmetry  evolves 
as 


9  E  +  4*2  E  =  9  B 

3c.  r  c  r  3C  9 


(68) 


If  avalanching  is  to  be  effective,  the  azimuthal  magnetic  field  changes 
little  during  the  breakdown  process.  Dropping  the  right-hand  side  of 
Eq .  (68)  and  setting  the  total  electric  field  E  =  E  allows  us  to  combine 
and  integrate  Eqs.  (67)  and  (68)  to  obtain 
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4n<T(r,t)  *  4*^(0  +  p  J  ^  v(E/p).  (69) 

E 

Here  o^( r)  -  <j(r,0)  is  the  preexisting  conductivity,  E(r,C)  is  the 
instantaneous  electric  field,  and  EQ(r,C)  is  the  effective  initial  field 
given  by 

Eo(r,t)  =  B0(r,O  =  2I(r,t)/rc.  (70) 

Because  BQ  is  not  truly  constant,  Eqs.  (69)  and  (70)  actually  represent 
only  an  approximate  upper  bound  for  the  conductivity  a  outside  the  beam. 

The  equations  are  fairly  accurate,  however,  in  the  region  of  greatest 
interest  where  4nar/c  «  1. 

If  the  initial  preionization  is  small,  the  final  conductivity  after 
breakdown  is  given  by 

E 

o 

4n<if  *  p  J  5(E/P)  »  4^  (71) 

Eb 

where  the  breakdown  field  strength  in  air  is  E^/p  e  30  kV/cm-atm.  That 
is  largely  independent  of  o\  indicates  that  preionization  alone  is  unlikely 
to  produce  a  well-defined  hollow  conducting  channel  suitable  for  tracking. 

Potentially  more  effective  is  a  low-density  channel  for  which  the 
product  pr  is  small  within  the  channel  but  large  everywhere  outside  the 
channel  and  beam.  If  v  is  a  strong  function  of  E/p,  as  usually  true,  a ^  is 
large  only  within  the  channel.  According  to  Eq.  (71),  the  channel 
conductivity  produced  by  avalanching  is  then  limited  to 


(72) 


4n ac  <  4nof(ac,C)  <  pc  v(Eq/pc) 

j  where  pc  is  the  channel  density  and 

Eq  .  Eo(ac,C)  *  2I(ac,C)/acc.  (73) 

To  obtain  high  conductivity  satisfying  condition  (12)  over  a  broad  channel 
thus  typically  requires 

I 

■ 

|  p  v(E  /p  )  »  c/a  .  (74) 

|  coc  c 

Condition  (74)  together  with  Eq.  (73)  defines  an  allowed  range  for  the 

channel  density  p  in  terms  of  the  avalanche  function  \),  the  inner  channel 

I  c 

radius  and  the  net  current  I(ac,C)  evaluated  at  the  pinch  point  Cp. 

I  The  maximum  channel  density  allowed  by  condition  (74)  depends  sensitively 

on  the  pinch-point  current  I(ac»Cp)  which  evolves  slowly  as  the  beam  head 
erodes.  The  ambient  gas  density  at  r  <  a£  should  satisfy 

I  Pn  >  2  P  (75) 

,  o  -  c 

'  to  produce  a  well-defined  channel. 

The  use  of  avalanching  to  produce  hollow  conducting  channels  has  a 
!  second  important  advantage  that  is  not  readily  apparent  from  the  preceding 

I 

!  axisymmetric  analysis.  A  displaced  beam  introduces  azimuthal  asymmetry 

into  the  channel  space-charge  fields.  This  asymmetry  is  amplified  by  the 
avalanche  function  v(E/p)  to  produce  a  channel  conductivity  a peaked 
'  nearest  the  beam.  The  asymmetry  in  a c  enhances  the  asymmetry  in  the 

i  channel  return  current  and  thus  enhances  the  channel  tracking  force.  The 

j  nonlinear  character  of  v  can  produce  a  strong  asymmetry  in  ctc  that  could 

i 

\  raise  the  tracking  force  by  more  than  a  factor  of  two.  At  late  times,  the 

t 
i 
» 


i 

[ 
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force  remains  tracking.  Even  channels  that  are  weakly  preionized  but  not 
rarefied  may  be  capable  of  producing  modest  tracking  forces  via  asymmetric 
avalanching.  Asymmetric  avalanching  is  likely,  however,  to  degrade  beam 
stability. 

Although  channel  rarefaction  can  produce  strong  tracking,  it  is  not  an 
energy-efficient  process.  The  hydrodynamic  energy  required  to  halve  the 
ambient  gas  density  is  orders  of  magnitude  higher  than  the  ionization 
energy  required  to  create  the  necessary  channel  conductivity  directly. 
Tracking  experiments  in  rarefied  channels  of  modest  length  may  nonetheless 
be  practical  if  direct  gas  ionization  techniques  are  not  readily 
available. ^ 


5.  Other  Effects 

We  describe  here  several  other  effects  that  can  influence  tracking. 
The  first  is  beam  rise  time.  If  the  beam  current  rises  slowly,  it 
generates  conductivity  gradually  so  that  the  pinch  point,  as  defined  by 
condition  (14),  occurs  late  in  the  beam  pulse.  If  the  pinch  point  occurs 
at  a  location  Cp  >>  CTC^*  the  tracking  force  is  reduced  by  the  fraction 
(cxcj/Cp)  <<  1,  as  suggested  by  Eqs.  (61)  and  (62).  To  minimize  this 
reduction,  the  channel  dipole  decay  time  should  therefore  satisfy 


ct  .  >  C  • 
cl  p 


(76a) 


This  imposes  an  additional  constraint  on  either  or  both  the  channel 
resistance  Rc  and  conductivity  ^co»  as  required  by  either  result  (37)  or 
(62).  For  a  pencil  beam  with  a  constant  rise  rate  i^,  the  pinch  point 
occurs  at  Cp  <  c  (a^/2kljj)  .  Condition  (76a)  is  thus  typically  satisfied 
provided 


Tcl  >  <abp/2kib> 


1/2 


(76b) 


Here  is  the  beam  radius  at  the  pinch  point,  and  k  is  the  direct- 
production  coefficient  defined  through  the  relationship 


3  a, 

ar  -  kJb  + 


(77) 


-4 

In  nitrogen,  k  =  5  x  10  cm/statcoulomb. 

Even  if  condition  (76a)  is  satisfied,  the  tracking  force  at  Cp  still 
falls  with  increasing  rise  time  because  the  beam  current  at  £  falls.  For 
example,  for  a  pencil  beam  with  constant  i^, 


Kac.Cp)  <  !„«„>  <  (ib  ab/2k)1/2. 


(78) 


The  decrease  in  the  tracking  force  is  partially  offset,  however,  by  an 
increase  in  the  tracking  duration.  Moreover,  a  short  beam  rise  time  (i.e., 
large  i^)  can  degrade  tracking  by  initiating  avalanching  in  the  inductively 
generated  axial  field  E^.  On-axis  avalanching  causes  afe  to  rise  rapidly  so 
that  the  pinch  point  moves  forward  and  I^CCp)  and  the  tracking  duration  are 
reduced.  Optimal  beam  rise  time  for  maximizing  the  tracking  force  at  the 
pinch  point  is  thus  expected  to  occur  when  the  on-axis  field  Ez(0,Cp)  -  E^ 
where  E^  is  the  gas  breakdown  field  strength. 

A  second  influence  on  tracking  is  a  dipole  chemistry  effect  arising 
from  field  dependence  in  the  plasma-electron  mobility  p.  Hubbard,  et  al. 
first  discovered  this  effect  in  hose  studies.  A  usually  adequate 
expression  for  the  conductivity  is 

a  =  en  u  (79) 

e 


where  ng  is  the  plasma-electron  density.  The  mobility  is  defined  by 
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(80) 


U  -  vd/E 


where  vd(E)  is  the  plasma-electron  drift  velocity.  Straightforward 

linearization  of  the  plasma  current  density,  J  -  oE  »  en  v.,  reveals  that 

*  p  e  a 

the  dipole  equations  (23)  and  (31)  are  more  accurately  represented  by 
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and 
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respectively,  where  the  effective  "dipole"  conductivities  are  defined  by 
abl  *  enebWd  <83> 

and 


cl 


en  u  . . 
ec  d 


(84) 


Here  nec  is  the  plasma-electron  density  in  the  channel,  n^  is  the  plasma- 
electron  density  generated  symmetrically  about  and  by  the  beam,  and  is 
the  differential  mobility  given  by 


wd 


1 


aln<y 

31n(E) 


W 


J 


(85) 


where  vm  »  eE/n^v^  is  an  average  momentum- transfer  collision  frequency  for 
the  plasma  electrons.  Not  included  in  these  equations  are  azimuthal 
asymmetries  in  n  ^  caused  by  avalanching  as  discussed  in  the  previous 
section. 

The  principal  change  introduced  by  Eqs.  (81)  and  (82)  is  that  the 
dipole  decay  times  T  ^  and  are  determined  not  by  the  conductivities 

and  <x,  but  by  the  effective  conductivities  a  ,  and  o,,.  If  the  collision 

b  1  cl  bl 

frequency  increases  with  field  E,  as  is  typical  of  most  weakly  ionized 

gases,  the  dipole  decay  times  are  reduced.  The  reduction  is  modest, 

however,  in  gases  such  as  nitrogen  where  0.5  <  w^/y  <  1.  Increasing  the 

channel  ionization  by  a  factor  of  two  or  so  is  thus  often  sufficient  to 

insure  that  condition  (15)  for  t  .  remains  satisfied.  In  unusual*^  or 

cl 

highly  ionized  gases,  y^  can  exceed  y  and  the  dipole  decay  times  are 
increased.  In  some  gases,  y^  can  be  negative  but  the  plasma  currents  are 
then  inherently  unstable.  In  all  cases  the  tracking  force  never  exceeds 
the  value  predicted  by  Eq.  (9). 

A  third  effect  on  tracking  is  channel  fill-in.  For  some  applications, 
channel  fill-in  cannot  be  avoided,  and  hence  a  truly  hollow  channel 
configuration  cannot  be  maintained.  Hollow  channels  produced  by  laser 
preionization,  for  example,  suffer  from  laser  refraction  and  diffraction 
which  cause  the  laser  beam  to  optically  distort. 

The  effect  of  channel  fill-in  on  tracking  can  be  assessed  by  adding  a 
constant  floor  conductivity,  o for  r  <  a^.  If  the  fill-in  conductivity 
is  low  (4nCT^ac/c  <  1  and  <<  aco),  channel  tracking  is  probably  improved 
because  weakens  or  eliminates  deleterious  electrostatic  effects.  If 
is  large,  however,  it  produces  an  overlap  detracking  force  given  according 
to  Eq.  (31)  by 
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(86) 
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is  an  average  axial  electric  field.  Because  the  monopole  field, 

Ezq  *  -  3Aq/3C,  usually  varies  slowly  with  r,  the  average  field  satisfies 

Ezo  *  [l(ac,0  -  Ib(C)]/[naJob  +  na*of].  (88) 

If  we  define  the  peak  channel  tracking  force  by 

Ft  =  2eftrI(ac,C)/acc?  (89) 

we  find  that 


Ff 

fyo  -  l(ac,0) 

f  . 1 

f<vt>  j 

2  2 
(*  ff f  ♦  a.ff,J 

(90) 


This  equation  indicates  that  channel  fill-in  does  not  appreciably  degrade 
tracking  unless  it  causes  the  net  current  inside  the  channel  to  fall  below 

I(ac,C)  <  Ib(C)/2.  (91) 


In  that  case,  the  net  force  can  become  detracking  and  eject  the  beam  from 
the  channel. 

At  beam  injection,  condition  (91)  is  satisfied  and  the  net  channel 
force  is  detracking,  provided  4najac/c  >  1.  The  plasma  currents  in  that 
cause  detracking  then  relax  while  the  net  current  I(ac»C)  grows.  For 
-*  0,  the  characteristic  relaxation  time  for  the  detracking  force  to 
decay  and  I(acC)  to  grow  is  given  by 


3! 


As  I(ac>C)  grows,  the  tracking  force  from  dipole  currents  in  a 
develops.  The  development  of  this  tracking  force  is  hampered,  however,  by 
its  relaxation  on  a  characteristic  time  x  j.  Only  if 


Tcl  >  Tf 


(93) 


can  the  tracking  force  fully  develop  and  overcome  the  detracking  force 
caused  by  plasma  currents  in  a According  to  Eqs.  (62)  and  (92),  a  broad 
hollow  channel  of  conductivity  acQ  can  therefore  produce  strong  tracking 
only  if  the  fill-in  conductivity  satisfies 


ac  <  a  /2. 
f  co 


(94) 


Conversely,  if  the  fill-in  conductivity  exceeds  fifty  percent  of  the  off- 
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axis  peak  conductivity,  detracking  is  likely.  Hubbard  and  Slinker  have 
performed  numerical  simulations  that  confirm  this  prediction. 

A  fourth  effect  that  can  degrade  tracking  is  azimuthal  nonuniformity 
in  the  channel.  Nonuniformity  caused  by  laser  hot  spots,  preferential 
avalanching,  or  other  means  should  not  be  too  disruptive  provided  the 
nonuniformity  is  mild.  The  principal  effect  is  to  displace  the  effective 
channel  center  and  to  produce  detracking  in  the  beam  body,  much  as  occurs 
in  the  magnetic  cage.  An  example  of  strong  nonuniformity  is  a  three- 
conductor  or  four-conductor  cage  with  one  conductor  missing.  In  that  case 
the  absence  of  conductivity  over  a  broad  region  would  lead  to  violent 
ejection  of  the  beam. 
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A  dynamical  effect  that  invariably  degrades  tracking  is  beam 

expansion.  Gas  scattering,  ohmic  losses,  and  finite  beam  emittance  all 

cause  the  beam  head  to  expand  and  flare  radially  outward  as  it  advances. 

This  flaring  pushes  the  pinch  point  Cp  backward  so  that  the  beam  head 

18  19 

continually  erodes.  ’  Because  the  channel  tracking  force  is  of  finite 
magnitude  and  extent,  any  given  beam  slice  experiences  a  positive  tracking 
force  for  only  a  finite  time  before  erosion  sweeps  past.  If  the  duration 
is  too  short  to  allow  the  beam  to  respond,  tracking  is  curtailed.  The 
curtailment  is  small  only  if  the  erosion  rate  satisfies 


d  Cn  AC 

t— ^  <  — - 

dz  zt 


(95) 


Here  characterizes  the  extent  of  the  tracking  region,  and 


1  1 
zt  =  j  (mc2)2[y/Ft(y)]2 


(96) 


characterizes  the  axial  distance  required  for  the  tracking  force  Ft  to 
deflect  the  beam.  For  hollow  channels  with  ACt  =  erosion  can 

typically  be  neglected  only  if 
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where  IA  =  0ymc  /e  is  the  Alfven  current.  Erosion  is  usually  most  severe 
for  lov-energy  beams  and  often  seriously  threatens  channel  tracking. 

A  related  problem  is  the  rise  in  the  beam  radius  a^  at  the  pinch 
point  as  the  beam  nears  the  end  of  its  range.  Unless  the  channel  radius  ac 
increases  with  propagation  distance  z,  a^  will  eventually  exceed  ac>  The 
beam  is  then  likely  to  be  ejected  from  the  channel  or  suffer  rapid  erosion. 
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Loss  of  the  beam  is  inevitable.  Such  loss  is  generally  of  little 
consequence,  however,  because  the  loss  occurs  only  after  the  beam  has  lost 
most  of  its  energy  and  has  begun  to  rapidly  expand. 

6 .  Numerical  Simulations 

20 

Preliminary  numerical  simulations  using  the  DYNASTY  field  solver 
have  partially  confirmed  the  preceding  analytical  predictions.  Typical 
results  for  the  net  tracking  force  are  presented  in  Fig.  3.  In  these 
simulations  the  beam  was  displaced  by  0.05  cm  from  the  channel  center  with 
a  beam  current  density  given  by 

Jb(r,C)  =  (Ibo/Jiab)exp(-r2/ab)tanh(C/Cr).  (98) 

In  all  runs  avalanching  was  turned  off,  to  aid  comparison  with  theory,  and 
the  beam  parameters  were:  peak  current  I.  =  10  kA,  radius  ab  =  0.5  cm, 
rise  length  =  15  cm,  and  direct-production  coefficient  k  =  8.8  x  10-4 
cm/statcoulomb. 

In  curves  a  and  b  of  Fig.  3,  the  preexisting  conductivity  for  the 
annular  channels  was  given  by 

°c(r)  =  acoexpl_4(r_rc)2/5rc1  ("> 

with  6rc  =  rc  =  1  cm.  The  use  of  steeper  profiles  for  a ■  led  to  numerical 
difficulty.  Despite  the  presence  of  both  substantial  channel  fill-in  and 
substantial  beam-channel  overlap,  the  tracking  force  was  large.  In  curve  a 
where  4nocoa^/c  =  5  (with  Rc  =  6  S/cm  and  ct^  =  2  cm),  a  peak  tracking 
force  of  160  Gauss  per  cm  of  beam  displacement  was  reached  at 
C  =  5  cm.  Detracking  occurred  beyond  C  =  18  cm,  but  this  should  be  offset 
by  strong  body-head  coupling.  Indeed,  channel  fill-in  caused  the  pinch 


point  Cp  ■*  0  in  this  example.  Raising  the  channel  conductivity  to 
411^^^/c  «  50  (with  Rc  «  0.6  Q/cm  and  ct^  a  20  cm)  caused  the  tracking 
force  in  curve  b  to  continue  rising  until  it  peaked  at  575  Gauss  per  cm  at 
C  =  17.5  cm. 

Curves  c  and  d  were  derived  for  a  magnetic  cage  consisting  of  four 
solid  channels,  each  displaced  symmetrically  1  cm  from  the  cage  axis  and 
each  with  a  Gaussian  radius  of  0.5  cm.  The  net  cage  resistance  in  curve  c 
was  Rc  «  6  Q/cm  while  in  curve  d,  Rc  a  0.6  Q/cm.  Comparison  of  curves  a 
and  c  and  curves  b  and  d  demonstrates  that  the  four-conductor  cage  degraded 
the  tracking  force  produced  by  an  equivalent  smooth  annulus  by  <  20%. 

In  curve  e,  a  low  channel  conductivity  of  4rta  a. /c  =  1  was  used, 

co  b 

spread  over  a  broad  annulus  from  r  =  0.5  to  5  cm.  The  tracking  force  was 
reduced  but  still  large  and  peaked  at  36  Gauss/cm.  The  tracking  force 
persisted  for  a  surprisingly  long  time.  Whether  this  persistence  is  an 
artifact  of  the  code  or  accurately  represents  the  diffusive  decay  of  the 
tracking  force  remains  unclear. 


7.  Summary 

The  preceding  analysis  and  simulations  have  shown  that  hollow  channels 
of  even  modest  conductivity  can  produce  strong  tracking  forces  on 
relativistic  electron  beams  propagating  inside.  Such  channels  have  the 
added  advantage  of  improving  stability  in  the  beam  head.  The  combined 
traits  of  strong  tracking  and  improved  beam  stability  make  hollow  channels 
particularly  attractive  for  both  long-range  and  short-range  guidance  of 
charged  particle  beams. 

The  principal  requirements  for  successful  hollow-channel  tracking  are 

high  conductivity  and  conductance  in  a  channel  that  lies  outside  the  beam. 

Dipole  chemistry  effects,  slow  beam  rise,  and  moderate  channel  fill-in  or 

nonuniformity  may  weaken  the  tracking  force  but  rarely  eliminate  it.  The 

principal  practical  obstacles  are  creating  the  channel  and  finding  a 

stable,  slowly  eroding  mode  of  beam  propagation. 

Although  preliminary  numerical  simulations  have  partially  confirmed 

the  analysis,  further  numerical  work  is  required.  Static  simulations  are 

needed  to  assess  the  impact  of  avalanching  on  rarefied  or  weakly  preionized 

channels,  and  to  study  the  effects  of  dipole  chemistry,  beam  rise,  and 

channel  fill-in  and  nonuniformity.  Dynamical  simulations,  such  as  those 
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begun  by  Hubbard  and  Slinker,  are  needed  to  study  beam  stability,  beam 

erosion,  loss  of  beam  particles  that  lie  outside  the  annulus,  and  the 

behavior  of  the  beam  once  the  pinch  point  expands  outside  the  channel. 

We  mention  in  closing  that  hollow-channel  tracking  has  been 

2 

successfully  demonstrated  in  the  laboratory  by  Leifeste,  et  al.  and  by 

14 

Raleigh  and  Fernsler.  Leifeste,  et  al.  used  a  hollow  Nd:YAG  laser  to 


preionize  an  annular  channel  in  nitrogen  doped  with  diethylanaline. 


Successful  guidance  and  beam  transport  were  achieved  at  a  gas  pressure  of 
1  Torr  with  only  moderate  beam  hose.  More  recently,  Raleigh  demonstrated 
beam  guidance  and  partial  stabilization  at  pressures  up  to  40  Torr  in  air 
using  a  cage  constructed  of  four  resistive  rods. 

Acknowledgements 

The  authors  thank  Dr.  Richard  Hubbard  and  Dr.  Steven  Slinker  for 
sharing  their  numerical  results  on  hollow-channel  tracking  and 
stabilization.  The  comments,  suggestions,  and  encouragement  offered  by  Dr 
Martin  Lampe  are  gratefully  acknowledged  as  well.  This  work  was  supported 
by  the  Defense  Advanced  Research  Projects  Agency  under  ARPA  Order  No.  4395 
Amendment  No.  63,  and  monitored  by  the  Naval  Surface  Veapons  Center. 


References 


R.  F.  Hubbard,  R.  F.  Fernsler,  S.  P.  Slinker,  A.  V.  Ali,  M.  Lampe,  G. 
Joyce  and  J.  M.  Picone,  "Stability  of  Relativistic  Electron  Beams  in 
Density  Channels  I",  5th  International  Conference' on  High-Power 
Particle  Beams,  San  Francisco,  CA  (Sept.,  1983). 


H.  C.  Chen  and  H.  S.  Uhm,  Bull.  Am.  Phys.  Soc.  29,  1275  (1984). 


2. 

G 

T 

Bull. 

3. 

H. 

C. 

4. 

P. 

A. 

P. 

Vr 

5. 

J. 

R. 

Vitkovitsky  and  A.  W.  Ali,  Phys.  Rev.  Lett.  41,  174  (1978);  D.  P. 
Murphy,  M.  Raleigh,  E.  Laikin,  J.  R.  Greig  and  R.  E.  Pechacek, 
"Electron  Beam  Transport  through  the  Atmosphere  in  Reduced-Density, 
Current-Carrying  Channels",  Proceedings  of  the  Ninth  International 
Symposium  on  Engineering  Problems  in  Fusion  Research,  Chicago, 
Illinois,  1981  (IEEE  Publication  81CH1715-2  NPS  Vol.  II,  p.  1548). 

J.  Benford  and  B.  Ecker,  Phys.  Rev.  Lett.  26,  1160  (1971). 

W.  E.  Martin,  G.  J.  Caporaso,  V.  M.  Favley,  D.  Prosnitz  and  A.  G. 

Cole,  Phys.  Rev.  Lett.  54,  685  (1985). 

Result  (2)  is  derived  by  applying  the  ultrarelativistic  frozen 
approximation  (19)  and  the  Lorentz  gauge  condition  to  Ampere's  law.  No 
further  approximations  or  assumptions  are  required. 

E.  P.  Lee,  "Calculation  of  a  Tracking  Force",  Lawrence  Livermore 
National  Laboratory,  UCID-19674  (1983)  [  DE83006669  ]  ;  R  J  Briggs,  private 
communication . 

Expression  (14)  is  used  for  convenience  to  characterize  when  magnetic 
forces  begin  to  overtake  electrostatic  forces.  Here  a.  is  the 
conductivity  evaluated  along  the  beam  axis.  Other  authors  define  the 
pinch  point  differently,  depending  upon  the  application. 

E.  P.  Lee,  Phys.  Fluids  21,  1327  (1978). 

M.  Lampe,  W.  Sharp,  R.  F.  Hubbard,  E.  P.  Lee  and  R.  J.  Briggs,  Phys. 
Fluids  27,  2921  (1984). 

R.  F.  Hubbard  and  S.  P.  Slinker,  private  communication. 

M.  Raleigh  and  R.  F.  Fernsler,  "Cage  Mode  Propagation  Experiments", 

NRL  Memorandum  Report  5784  (1987). 


D.  P.  Murphy,  M.  Raleigh,  R.  E.  Pechacek  and  J.  R.  Greig,  "Stability 
of  Relativistic  Electron  Beams  in  Density  Channels  II",  5th 
International  Conference  on  High-Power  Particle  Beams,  San  Francisco. 
CA  (Sept.,  1983). 

R.  F.  Hubbard,  S.  P.  Slinker,  M.  Lampe,  G.  Joyce  and  R.  F.  Fernsler, 
Bull.  Amer.  Phys.  Soc.  29,  1197  (1984).  See  also  R.  Hubbard,  et  al.t 
"Beam  Propagation  Studies  at  NRL,  Vol.  2",  NRL  Memorandum  Report  5412 
(1984).  ADA148198 

T.  G.  Christophorou,  S.  R.  Hunter,  J.  G.  Carter  and  R.  A.  Mathis, 
Appl.  Phys.  Lett.  41,  147  (1982). 

E.  P.  Lee,  "Model  of  Beam  Head  Erosion",  Lawrence  Livermore  National 
Laboratory,  UDIC-18768  (Aug.  8,  1980). 

W.  M.  Sharp  and  M.  Lampe,  Phys.  Fluids  23,  2383  (1980). 

B.  Hui  and  M.  Lampe,  J.  Comp.  Phys.  55,  328  (1984). 


Figure  1.  Hollow-Channel  Tracking.  A  displaced  electron  beam  propagates 

in  an  annular  conducting  channel  bounded  by  the  heavy  concentric 
circles.  A  symmetric  ring  of  beam  current  (dashed  circle  in  a) 
entirely  inside  the  channel  experiences  a  net  tracking  force 
initially  and  zero  force  later.  A  beam  ring  (b)  entirely 
outside  the  channel  always  experiences  zero  net  force.  A  beam 
ring  (c)  intersecting  the  channel  always  experiences  a  net 
detracking  force  caused  by  channel  return  current  flowing 
outside  the  ring  (shaded  area). 


Figure  2.  Four-Conductor  Magnetic  Cage.  Four  solid  conducting  (shaded) 
channels  are  arranged  in  a  square  configuration  to  guide  an 
electron  beam  (dashed  circle). 
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Figure  3. 
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